Cosmological constraints on a decomposed Chaplygin gas 
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Any unified dark matter cosmology can be decomposed into dark matter interacting with vacuum 
energy, without introducing any additional degrees of freedom. We present observational constraints 
on an interacting vacuum plus dark energy corresponding to a generalised Chaplygin gas cosmology. 
We consider two distinct models for the interaction leading to either a barotropic equation of state 
or dark matter that follows geodesies, corresponding to a rest-frame sound speed equal to the 
adiabatic sound speed or zero sound speed, respectively. For the barotropic model, the most stringent 
constraint on a comes from the combination of CMB+SNIa+LSS(m) gives a < 1.44 x 10 at the 
95% confidence level, which indicates that the barotropic model must be extremely close to the 
ACDM cosmology. For the case where the dark matter follows geodesies, perturbations have zero 
^v^j ' sound speed, and CMB+SNIa+glSW then gives the much weaker constraint —0.1522 < a < 0.2556 

at the 95% confidence level. 

I. INTRODUCTION 

The dominance of dark energy and dark matter in the present cosmological energy density has led many 
/"*\ authors to seek a unified description in terms of a single dark component, referred to as unified dark matter or 
quartessence 0], which should explain both the current accelerated expansion of the universe and the role of 
^■•^ . non-baryonic dark matter in structure formation. One candidate from this kind of model, which has attracted 
' a l°t °f attention, is the generalized Chaplygin gas (GCG) model Q. It is a fluid, defined by an exotic equation 
Q^i, of state: 

O '. A 

£■ P gCg = --^, (1) 

ryj ' PgCg 

a ■ 

i 1 with model parameters a and A. The GCG model with a = 1 reduces to the original Chaplygin gas model 

0. Predictions of the GCG cosmology has been tested against the different observational data, including 
type la supernovae 0-0, * ne positions of cosmic microwave background radiation peaks 0-0], the X-ray gas 
mass fraction of clusters 0], Hubble parameter- redshift data [ill E3; the lookback time-redshift data fl3j . 

■ gamma-ray bursts [3] , and as well as various combinations of the updated data 0] . 

(f") | Going beyond constraints provided by the homogeneous background cosmology, the inclusion of inhomo- 

■ geneous linear perturbations in the GCG has an effect on both the cosmic microwave background (CMB) 
anisotropics and structure formation. In particular, the imprint of a large integrated Sachs- Wolfe (ISW) signal 

16], which arises from a time-dependent gravitational potential when dark energ y d ominates at late times in 
the Universe, could be observed on the CMB anisotropies 0, 0] ■ Recently, in [191 ] two of us revisited con- 
straints on the GCG model from the combination of CMB, baryon acoustic oscillations (BAO), and type la 
supernovae (SNIa), including the full CMB temperature and polarization power spectra data. A tight bound 
\ for the parameter a was obtained up to the order of 10 -3 . 

The effect of the GCG on structure formation also contains important information. In [20| . it was argued 
that in order to predict large-scale structure (LSS) consistent with the data from the 2dF survey 0], the GCG 
model has to be indistinguishable from ACDM, with a < 0(1O~ 5 ). The slightly looser constraint, a < O(10~ 4 ) 
was obtained in 0] by fitting the predicted baryon matter power spectrum, rather than the total matter power 
spectrum, to the Sloan Digital Sky Survey (SDSS) DR7 power spectrum data. In either case a stringent limit 
on the a parameter is required from the LSS power spectrum data because there are instabilities or oscillations 
of GCG power spectra due to the adiabatic pressure perturbation produced for non-zero a. This problem can 
be avoided in some unified dark matter models by requiring a fast transition between a CDM-like era and a 
ACDM-like phase 0,113. 

In order to save the GCG model some type of intrinsic entropy perturbation must be introduced in the GCG 
fluid to counteract the effect caused by the adiabatic GCG pressure perturbation, as in the so-called silent 
GCG model [25[. The combined constraints from data on both the silent GCG and standard GCG model are 
presented in [26| . It is shown that with a vanishing GCG pressure perturbation there is a much wider allowed 
parameter range, —0.3 < a < 0.7 0]. Subsequently, the ISW data was used to constrain the silent GCG 
model 0. The required entropy perturbation can be introduced by decomposing the GCG fluid into two 
components of dark matter and dark energy [29l |. but the possible ways of decomposition are not unique. In 
[30l ] the GCG model was split into dark matter and vacuum energy, considering a homogeneous vacuum when 
studying density perturbations in the Newtonian limit. 
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Recently it has been shown that any unified dark matter model can be decomposed into dark matter inter- 
acting with a vacuum energy (with equation of state P = —p) without introducing any additional degrees of 
freedom [3l| . In general an interacting vacuum ener gy is inhomogenous, and a description of the dynamics of 
inhomogeneous vacuum energy was presented in [3lLl32 |, 

In this paper we will study in detail the constraints on a decomposed GCG model described by dark matter 
interacting with inhomogeneous vacuum energy. A closed set of perturbation equations are obtained once a 
covariant form is specified for the interaction. Firstly, we consider a model in which the vacuum energy is 
a function of the local matter density. We show this leads to a barotropic equation of state for the dark 
mat ter+ vacuum, and hence reproduces previous results for a unified GCG. Secondly we consider an energy- 
momentum transfer parallel to the matter four-velocity, leading to a geodesic flow for the matter. In both cases 
we consider background cosmologies that mimic the GCG cosmologies, but the evolution of linear perturbations 
differs due to the presence of non-adiabatic pressure perturbations. We pay particular attention to their imprints 
on the CMB and LSS power spectra. To constrain the barotropic and geodesic model parameters we will use 
different combinations of current cosmological data including the CMB temperature and polarization power 
spectra (33[, SNIa from the Union2.1 compilation of the Supernova Cosmolo gy P roject Collaboration [34fl . the 
BAO distance measurements (35[, the LSS power spectrum from SDSS DR7 |36j. and the galaxy-ISW (gISW) 
cross-correlation power spectra from [37j | . 

This paper is organized as follows. In the next section, we review the background equations for the GCG 
model. In Sec. Ill, the linear perturbation equations for both the unified GCG model and the two decomposed 
models are presented. In Sec. IV, we show the effects on CMB and LSS power spectra in the two decomposed 
models, then perform a global fitting to observational data by using the MCMC method and discuss the resulting 
constraints on the interaction. In the final section we present our conclusions. In what follows (unless otherwise 
specified) the illustrative plots are shown using the joint mean values for cosmological parameters from WMAP7 
data [H: n b h 2 = 0.02255, tt dm h 2 = 0.1126, h = 0.702, r = 0.088, n s = 0.968, and A s = 2.43 x 10" 9 at k s0 = 
0.002MPC" 1 . 



II. CHAPLYGIN GAS COSMOLOGY 



We will consider a spatially-flat, homogeneous and isotropic Friedmann-Robertson- Walker universe with met- 
ric 

ds 2 = -dt 2 + a 2 (i)%cfeW , (2) 
where the Friedmann equation is given by 
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H 2 = \^-J =—(Pb + Pr + P gCg ), (3) 

where p b , p r and p g Cg are the energy densities of baryons, radiation and a generalised Chaplygin gas (GCG) 

which can act as both a dark energy component, accelerating the universe at late times, and dark matter. The 
energy conservation equations for the various components read 

p b + 3Hp b = , (4) 

p r + 4H Pr = , (5) 

p gC g + 3if (p gC g + PgCg) = . (6) 

The equation of state for GCG H} then allows us to integrate Eq. ^ to obtain 



PgCg = PgCgO 



, n l/(l+a) 

B s + (l-B s )a~^ 1+a A , (7) 



where PgCgO is present value of energy density of GCG (when a = clq = 1) and we have replaced A in Eq. (TT|) 
with the dimensionless model parameter, B s = A/p g ^ . 

In a previous work we have shown that a unified dark matter fluid with an arbitrary equation of state can 
be decomposed into pressureless dark matter interacting with a vacuum energy, V, with an appropriate energy 
transfer [31] . For the GCG solution ([7]), we can write 

PgCg = Pdm + V , (8) 

such that, from Eq. ^ 

A = V( Pdm + V) a . (9) 



3 



The conservation equations of dark matter and vacuum energy are 

Pdrn + SHpdm = ~Q, 

V = Q. 

Thus from Eqs. ©, CGI]) and (HI]), we obtain 

PdmV 



Q = 3aH- 



Pdn 



V 



(10) 

(11) 



(12) 



In the decomposed model ([5]), the parameter, B s in Eq. ([7]), can be written in terms of the present values of 
the energy densities, B s = [V/(pd m + V)]o- The decomposed model is characterised by a single dimensionless 
parameter a [3li | , but an FRW solution is specified by two boundary conditions (such as the present density of 
matter and vacuum) whereas the original GCG model is specified by two model parameters, a and A (or B s ) 
and one boundary condition. This already suggests that the decomposed model may naturally accommodate 
non-adiabatic density perturbations, in addition to the adiabatic perturbations of the original GCG. As we shall 
see, this can have important consequences for the evolution of perturbations in a decomposed model. 

Note that the Friedmann equation ([3]) can be written as 



H 2 = H 2 



f2f,a + Q r a + (£V + ^dra) 



n 



din 



-3(l+a) 



Jly + fldm Qv + Qdm 

where the dimensionless fractional energy density of individual components is defined as 



(13) 



(14) 



which satisfies fit, + Q r + Qy + Qdm = 1- 

We illustrate the background evolution as a function of redshift for different values of the parameter a in 
Fig. [T] In the particular case when a = 0, the background evolution is identical to that for the ACDM model. 
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FIG. 1: The Hubble parameter as a function of redshift z for different values of parameter a, where the black line with 
a = corresponds to ACDM model. 



III. THE LINEAR PERTURBATION EQUATIONS 



We begin with a general scalar mode of the metric perturbation in an arbitrary gauge with the perturbed 
line of element [38H4ll 



ds 2 



The Einstein equations are [38H41I 



(1 + 2<j))dt z + 2ad l Bdtdx l + a 2 [(1 - 2i/>)5 ij + 2d l d j E] dx l dx 3 



3H (ip + H<j>j - —[ip + Ha) = -AirGSp , 

4, + H<j) = -4ttG(p + P)6, 

.... . / 2 V 2 

i> + 3HtP + H<j>+ (3H 2 + 2HU = AirG SP + -—=-11 

V 3 a z 



(15) 



(16) 

(17) 
(18) 

(19) 
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where a = a 2 E — aB is the shear, 5p is the total density perturbation, SP is the total pressure perturbation, II 
is the total anisotropic stress, and is the total covariant velocity perturbation . 

For each individual component, the covariant conservation equation allowing for energy-momentum transfer 
gives 

V,T[- = Qfo . (20) 

General relativity requires conservation of the total energy-momentum, V^T^ = 0, so we have = ^- 

The perturbed energy-momentum transfer four- vector is conventionally split into the energy transfer, Qj +SQi, 
and momentum transfer, //, relative to the total four- velocity [40l - |43j 

Q { £> = hQ/(l + 4>) - SQiAifi + Qj9)] ■ (21) 
The energy and momentum conservation equations for each fluid component are then (39l - l4ll [43[ : 

V 2 

5 Pl + 3H{S PI + SPj) -3(p/+P/)^ + (pi+ Pi) — (0i + o) = SQj + Qj^, (22) 



a* 

2 V 2 
3~ 2 



{pi + Pi)0i - 2>c 2 sI H( Pl + Pj)0! + ( PI + + SPj + = fi + QiO - (1 + cl^Qjh , (23) 



where c 2 sI — Pi /pi is the adiabatic sound speed. 

The numerical code CAMB [44j is written in the synchronous gauge, so in order toplace quantitative con- 
straints on the model parameters we will solve the equations in a synchronous gauge [451 ] where 

<t> = B = 0, ip = T), and k 2 E = -- - 3rj . (24) 

There is a residual spatial gauge freedom in the synchronous gauge which allows one to set the initial velocity 
of any component to zero. We see from Eq. (|23p that it is possible to work in a gauge in which 9j — at all 
times for non-interacting, pressureless matter. A comoving synchronous gauge is commonly used, in which the 
dark matter velocity is zero, 9d m = 0. This is no longer possible in models where the dark matter experiences 
a momentum transfer, fd m + Qdm{9 ~ Qdm) ^ 0. We can still work in a synchronous gauge, but note that the 
dark matter velocity is not necessarily zero. 

The time-time and time-space components of the perturbed Einstein field equations, (|16p and (|1T[) . in the 
synchronous gauge are 

5-5(£ ,+4 " a *)' (25) 

r) = -4vrG(p + P)6 . (26) 

Together with the energy and momentum conservation equations for each component, these form a closed set 
of equations. 

After decoupling, the baryon component is conserved independently and the perturbation equations for the 
baryon density contrast and velocity are given by 

b = 0. (28) 

Combining Eqs. (TT5|) . (fH?)) in the synchronous gauge and Eq. (|2"5")) . we have 

h + 2Hh = -8nG(5p + 35 P) . (29) 

By differentiating Eq. ((2"?]) with respect to time and combining Eq. (|28p and the above equation, we find a 
second-order differential equation for the baryonic density contrast 

5 b + 2H5 b = AnG(6p + 35P) . (30) 



1 Note that 9 = a(v + B) in terms of the total velocity potential v defined in, e.g., Ref. |4l| . 
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A. Unified GCG model 

In the unified GCG model, the energy-momentum tensor of the GCG is conserved and hence Eqs. (|22l) and 
give the perturbed energy and momentum conservation equations 

( k 2 h\ 

SgCg+3H(cl gCg -W g Cg)6 g Cg + (l+W g Cg)[--^0 gC g + -\ = 0, (31) 



gCg ~ 3c 2 s gCg H6 gC g + 5g C g - , (32) 



r 2 

1 + WgCg 

where the adiabatic sound speed for the GCG is fixed by the equation of state of the GCG 



c l,gCg = - aw a cg, (33) 

WgCg = -B a + {l-B a )a-XW (34) 

We can derive a second-order differential equation for GCG overdensity by differentiating Eq. (|3Tj) with 
respect to time and combining Eqs. (|2U1) and ([3"2"]) : 

SgCg + 2H[l+3(c 2 sgCg -w gC g)}8gCg 

+ [SH^gCg ~ WgCg) + 3H(ci gCg - WgCg)' + 9H 2 (cl gCg - Wg C gf + QH 2 {c 2 gCg - Wg C g)]Sg C g 

k 2 k 2 
= - c lgCg^ S gCg + (l + WgCg)3c 2 gCg H—e gC g + (l + WgCgHTrG{5p + 36P). (35) 

The first term on the right-hand-side of the above equation dominates in the small-scale limit and the adiabatic 
pressure perturbation causes the instabilities (or oscillations) of density perturbations for negative (or positive) 
a 0. 

B. Interacting vacuum energy + dark matter models 

In an interacting vacuum energy + dark matter model, the energy continuity equations (|22j) for the dark 
matter and vacuum energy in the synchronous gauge become 

( k 2 h\ 

fipdm + 3H5p dm + p drn I --^9 d rn + ^ J = &Qdm , (36) 

6V = SQ V , (37) 

while the momentum conservation equations (|23p are given by 

PdmQdm = j 'dm + Q 'dm(& ~ 6 'dm) , (38) 

-SV = fv + QvO. (39) 

Note that the vacuum energy has vanishing momentum and therefore the force exerted on the vacuum, di(fv + 
Qv@), is exactly cancelled by the gradient of the vacuum pressure, di(—SV). 

The conservation of total energy-momentum for the interacting dark matter and vacuum energy implies that 

Qv — —Qdm — Q ! $Qv — ~5Qdm , fv = —fdm ■ (40) 

We can then eliminate SQdm and fdm using the energy and momentum conservation equations for the vacuum 
(|37p and (|39p to leave the energy-momentum conservation equations for the dark matter 



•'./... " &dm 2® dm + o _ _ ~ (41) 



k 2 h SV 

—o"dm + 7T — 

a 2 p dn 

PdmOdrn = SV + QOdm, (42) 

where the background energy transfer Q is given by Eq. (|12p . We now need to specify the form of the vacuum 
interaction to determine the vacuum perturbation, SV, and obtain a closed set of equations. 
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1. Barotropic model 



Firstly, we consider a model in which the vacuum energy density is a function of local dark matter density 

V = V{p dm ) . (43) 

This implies that the energy transfer between the vacuum and the dark matter is along the gradient of the local 
dark matter density: Q v oc V 11 p dm . 
At linear order, we have 



SV 



V 



Sp. 



dm ■ 



(44) 



which shows that in this case the relative perturbations between the vacuum and the dark matter are adiabatic. 
The total pressure perturbation for the interacting dark matter+vacuum energy is 



SV = 



aV 



Pdrn + V 
2 



(Spdm + SV) 



,gCg( 6 Pdm + SV) . 



(45) 



Therefore, the sound speed for the interacting dark matter+vacuum energy is exactly the same as the adiabatic 
sound speed for the original generalised Chaplygin gas, c 2 gCg given by Eq. (|3"3"j) [HI, HI]. We show the effect on 
the sound speed caused by the parameter a in Fig. [2] 



FIG. 2: The sound speed as a function of scale factor a for the barotropic model 
In the synchronous gauge, we have from Eqs. (|4"Tj) . (l4"2l and (|4"4"|) 

Sdm 
8dm 



(l + a)Q , Pdm + (1 + a)V k 2 n _p dm + {l + a)Vh 

Odm i , t r 9 "dm 



Pdm + (1 + a)V 
Q „ aV 



Pdr, 



'dm 



p dm + (1 + a)V 



Pdm ~ 
Sdm ■ 



V 



Pdr, 



V 



2 ' 



(46) 
(47) 



It is clear that the barotropic model with a = reduces to the ACDM model at the linear order. Comparing 
with dark matter perturbations in the standard ACDM model, we can see that density contrasts and velocities 
of the interacting dark matter evolve very differently. In this case, due to the presence of momentum transfer, 
9dm will evolve to be nonzero, even if its initial value is set to zero. Therefore, the synchronous gauge is not 
comoving with the interacting dark matter and we calculate the perturbed equations in an arbitrary synchronous 
gauge [47j . 

Combining Eqs. (|46| and (j47]) with Eq. (|29p. we obtain a second-order differential equation for the dark 
matter overdensity 



-'dm 



2(1 + a)Q 



Pdm + (1 + a)V 

{l + a)Q 
p dm + (1 + a)V 



- 2H 



dm 



(l + a)Q 



+ {l + a)V 



2H 



(l + a)Q 



Sdn 



1 



aV 



Pdm + V J p dm a 2 



p dm + {l + a)V 
aV 



Idrn 



Pdr, 



V 



AirG(Sp + 3SP) . 



(48) 
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The first term on the right-hand-side of Eq. (|4"5]l dominates in the small-scale limit. Negative a causes insta- 
bilities of density perturbations. Conversely, there are oscillations of perturbations for positive a. In order to 
avoid the presence of negative sound speed, we will require a > in this case. 

We find that the perturbations for the barotropic model, Eqs. (l46l) . (l47l) and (|48|) . obey the same dynamical 
evolutions as those for the unified GCG model, Eqs. (I5T1) . (|52"j) and where we identify 



°gCg 



Pdn 



Pdm + (1 + ot)V 



y gCg 



(49) 



We should not be surprised that this interacting model reproduces the same behaviour as the GCG since by 
requiring V — V(pdm) we have set the pressure of the interacting matter+vacuum Pdm+v = —V to be an 
implicit function of pdm + V, i.e., we have a barotropic fluid. 

The evolution of the dark matter density perturbations for the barotropic model on a fixed scale is shown 
in Fig. [3l On the left panel, it is shown that increasing a causes the suppression of the growth of density 
fluctuations on the scale, k = 0.001[h Mpc -1 ]. For the smaller scale, k = 0.01[h Mpc -1 ], on the right panel, we 
can see that density fluctuations rapidly decay and oscillate with increasing a. 




FIG. 3: The dark matter density perturbations on fixed scales: k = 0.001 [h Mpc x ] (left panel) and k = 0.01 [h Mpc *] 
(right panel) as a function of scale factor for the barotropic model. The blue line, red line, green line, pink line and black 
line correspond to a = 0.1, 0.05, 0.01, 0.001, (ACDM model), respectively. 



2. Geodesic model 



An alternative model for the energy transfer is to consider an energy flow that is parallel to the 4- velocity of 
the dark matter: 

Q?dm) = -Qtf dm y (50) 

For this form of the covariant interaction there is no momentum transfer in the rest frame of dark matter 
[H, El|, I dm — Q(0 — 8 dm) = in Eq. (j2"Tj) . The velocity perturbation for dark matter is not affected by the 
interaction and obeys the standard equation, 0d m = from Eq. (|38[) . This allows us to work in a synchronous 
gauge that is comoving with the dark matter (8dm = 0). The dark matter 4- velocity, it M is thus a geodesic flow: 
M (dm)^A tM (dm) = ^- Hence we will refer to this form of energy-transfer as a geodesic model. 

The vacuum perturbation in the dark matter-comoving frame is identically zero, from Eq. (|42[) , 

6V com =8V + V8 dm = . (51) 

The pressure perturbation for the interacting vacuum+dark matter is therefore zero in the dark matter rest 
frame and thus the speed of sound for the interacting vacuum+dark matter is zero [H], H3 • Note that this 
implies that in this case there is a non-adiabatic pressure perturbation, due to the relative entropy perturbation 
between the dark matter and the vacuum energy [3lL H3| , in contrast to the barotropic model. 

From Eq. (|41[) . the perturbation equation for the dark matter density contrast in the rest frame of dark matter 
obeys 

5dm = ~ + —^—0~dm ■ (52) 
^ Pdm 
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We see that, again, the case with a — is completely equivalent to ACDM at linear order. 

A second-order differential equation for the dark matter overdensity can be derived from the above equation 
and Eq. (J2U), which gives 



'dm 



Pdr, 



2H S dm - 



Q 

2H— + 

Pdrn 



Q_ 

Pdm 



6 dm = 4ttG(S P + 36P) . 



(53) 



We plot the evolution of the dark matter density perturbations on different scales in Fig. 2] Compared with the 
growth in the ACDM model, positive a indicates that there is an energy transfer from dark matter to vacuum 
energy, so the growth of density perturbations with positive a is suppressed at early times during matter- 
domination. At late times, vacuum energy dominates in the universe, and the vacuum interaction drives the 
growth of dark matter density perturbations. Conversely, for the negative a the growth of density perturbations 
is enhanced at early times and is gradually suppressed at late times. Note that there is no instability associated 
with imaginary speed of sound for a < in the geodesic model. 




0.2 0.4 0.6 0. 
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FIG. 4: The dark matter density perturbations on fixed scales: k = 0.001 [h Mpc -1 ] (left panel) and k = 0.01 [h Mpc" 1 ] 
(right panel) as a function of scale factor for the geodesic model. The blue line, red line, black line, green line and pink 
line correspond to a — 0.2, 0.1, (ACDM model), -0.1, -0.2, respectively. 



IV. COSMOLOGICAL CONSTRAINTS 



A. Theoretical predictions of CMB and LSS power spectra 



First, we illustrate the imprint on the CMB and LSS power spectra for different values of the parameter a 
for both the barotropic model and geodesic model. For the barotropic model we will restrict the value of a to 
be greater than or equal to zero in order to avoid instabilities due to an imaginary sound speed, but we will 
extend a to negative values for the geodesic model since the sound speed remains real in that case. 

The effects on the CMB power spectra are shown in Fig. [5] From the left panel, positive a leads to an 
energy transfer from dark matter to vacuum energy. Increasing a causes the matter-radiation equality to 
happen earlier and hence suppresses the first acoustic peaks. Conversely, negative a delays the epoch of matter- 
radiation equality, and hence amplifies the acoustic peaks, as shown in the geodesic case on the right panel of 
Fig. [5] 

One source of CMB anisotropics at low multipoles is the ISW effect. The amplitude of the ISW signal on large 
scales is affected not only by the epoch of dark matter-vacuum energy equality, but is also modified directly 
by the evolution of the gravitational potential, which differs in the barotropic and geodesic models. These two 
competing mechanisms give the final ISW signal. 

In order to clearly express the ISW effect caused by the time- varying gravitational potential, we calculate the 
Newtonian potential, ip N , corresponding to the metric perturbation, ip, in the Newtonian or longitudinal gauge 

(B = E = 0) d3. 

This is given by the relativistic Poisson equation [|(| 

k 2 ip N = -AwGa 2 pA , (54) 
where the total comoving density perturbation is given by 



(55) 
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including the vacuum perturbation. The traceless, space-space component of the Einstein field equations, (|18l) . 
gives the relation between the two scalar potentials in the Newtonian gauge [4j| : 

■4, N - (j) N = 87rGn , (56) 

where II is related to total the anisotropic stresses. 

Now the ISW temperature anisotropy is given by the following source: 

s/sw = <y + i N 

= - 8 " G s(^ a+n )- < 57 > 

At late times, the relativistic particles (photons and neutrinos) can be neglected and there is no anisotropic 
stress. Therefore, it is the comoving density perturbation that drives the late-time ISW effect. Fig. [5] shows 
how the parameter a affects the time evolution of the gravitational potential in the two models. We see that for 
the same value of a we find a much larger variation in the Newtonian potential at late times in the barotropic 
model due to the non-zero sound speed, as shown in Fig. [3J This will be reflected in much tighter bounds on 
the value of a in the barotropic model coming from the overall CMB anisotropics, and these constraints become 
even tighter when the galaxy clustering data is also used to identify the late-time ISW effect. 

Note that the overall CMB anisotropies, including the ISW effect, is the same in the barotropic model and 
the original unified Chaplygin gas model since the evolution of the total matter+vacuum density perturbations 
is the same as the density perturbations in the unified GCG model. 




i i 

FIG. 5: CMB temperature power spectra vs multipole moment 1. The black solid line is for the ACDM model. Left 
panel: Other solid lines and dashed lines for the same values of parameter a correspond to the barotropic model and 
geodesic model, respectively. Right panel: The dashed line and dotted line correspond to a = 0.1 and a = —0.1 in the 
geodesic model, respectively. 

In Fig. we present the dark matter power spectrum, the baryonic power spectrum, and the total matter 
power spectrum (baryon and dark matter) for the barotropic model. The figure shows that, as expected, the 
dark matter has an oscillating power spectrum in the small-scale limit for non-zero a. The larger the parameter 
a, the larger the scale where the oscillations of dark matter power spectra begin. By comparison, the baryonic 
power spectra behave much more smoothly as they are only gravitationally coupled to the oscillating dark 
matter density [22| . The total matter power spectra combining both dark matter and baryon densities therefore 
shows intermediate behaviour, but cannot avoid the rapid decay of power spectrum on the small scales for 
non-zero a. 

Fig. [5] shows the total matter power spectra for the geodesic model. The problematic oscillations of the dark 
matter power spectra do not occur in this case since the sound speed is zero. The only oscillations seen are the 
baryon acoustic oscillations, also seen in ACDM. The matter power spectrum is enhanced on small scales for 
a > due to the earlier matter-radiation equality which moves the turnover in the matter power spectrum to 
smaller scales. 



B. Constraint method and data 

Next, we will use a variety of cosmological data to constrain these two models for the matter+vacuum inter- 
action, according to the Markov Chain Monte Carlo (MCMC) method, using the publicly available CosmoMC 
package [5(3, [EH- The cosmological parameter set is P = {rtbh 2 , fldmh 2 , 9s, t, a, n s , \og[10 10 A s }} , where fif,/i 2 
and Qdmh 2 are the physical baryon and dark matter densities, ®s is the ratio (multiplied by 100) of the sound 
horizon and angular diameter distance, r is the optical depth, a is the dimensionless generalised Chaplygin 
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FIG. 6: The evolution of the gravitational potential as a function of scale factor on the large scale k = 0.001[h Mpc -1 ]. 
The black solid line is for the ACDM model. Other solid lines are for the barotropic model, where the blue line, red 
line and green line correspond to a = 0.1, 0.2, 0.3, respectively. The dashed lines are for the geodesic model, where the 
blue line, red line and green line correspond to a — 0.1, 0.2, 0.3, respectively. The dotted line is the evolution of the 
gravitational potential for a — —0.1 in the geodesic model. 
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FIG. 7: The LSS power spectra as a function of wavenumber for the barotropic model. On each panel, the black line 
is the total matter power spectrum for the ACDM model. The pink line, red line and blue line are the baryonic power 
spectrum, the total matter power spectrum and the dark matter power spectrum, respectively. The separate panels (left 
to right and top to bottom) correspond toa= 10 -5 , 10 -4 , 10 -3 , 10 -2 , respectively. 
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FIG. 8: The total matter power spectra as a function of wavenumber for the geodesic model. The blue line, red line, 
black line, green line and pink line correspond to a = 0.2, 0.1, (ACDM model), -0.1, -0.2, respectively. 



gas model parameter, n s is the scalar spectral index, and A s is defined as the amplitude of the initial power 
spectrum. Note that the dimensional parameter A in the original definition of the Chaplygin gas (|T|) is given in 
terms of fldm and f2y using Eq. ([9]) as 



A — fly {Q<lm + ^V)' 



\8ttGJ 



l+a 



(58) 



and Qy = 1 — fife — fi r — £ldm, by the requirement of spatial flatness. 

We will use different combinations of data (CMB+SNIa+BAO, CMB+SNIa+LSS, and CMB+SNIa+glSW). 
For the LSS power spectra we will only use the linear matter power spectrum up to k ~ 0.1[h Mpc^ 1 ]. In 
addition, we use a top- hat prior for the cosmic age, i.e., lOGyr < to < 20Gyr, impose a weak Gaussian prior on 
the physical baryon density fi^/i 2 = 0.022 ± 0.002 from big bang nucleosynthesis [HI], and use the value of the 
Hubble constant Hq = 73.8 ± 2.4km s^ 1 Mpc _1 [53[ taking a Gaussian likelihood function. 

The theoretical galaxy-ISW (gISW) cross-correlated power spectrum, -Pj^,, where $ = ip N + cb N , can be 
calculated by a modified CAMB code. We modify the code by referring to ISiTGR [54], so that a power 
spectrum can be directly calculated from the individual transfer function. In the absence of magnification bias, 
the gISW cross-correlation is written as 



r s,T = T C 'MB 
1 (/ + 1/2) 2 



dzb x (z)U x (z)P s £ 



1 + 1/2 



(59) 



where the galaxy samples in [33] are collected from the Two Micron All Sky Survey (2MASS) Extended Source 
Catalog (XSC), which is divided into 4 bins: 2MASS0-3, the luminous red galaxies (LRGs) from SDSS split 
into 2 sample bins: LRG0-1, the photometric quasars from SDSS with 2 redshift bins: QSO0-1, and NRAO 
VLA Sky Survey (NVSS): NVSS. Namely, there are nine galaxy sample bins x. The scale-independent bias 
factor b x (z) is assumed to relate the observed projected galaxy overdensity to the matter density. H x (z) is the 
redshift distribution. x( z ) is the comoving radial distance. In ISWWLL code, the function f x (z) = b x (z)H x (z) 
is determined by fitting auto power spectra and cross power spectra between the samples: 



dzf x (z)f y (z)D 2 (z)P St s 



1 + 1/2 



(60) 



even if the methods differ for the various samples. 

For the 2MASS galaxies, the redshifts are identified by matching 2MASS galaxies with SDSS main galaxies. 
The nonlinear Q model [55j (denoted by the bold text to be distinguished from the background interaction Q) 
of the galaxy power spectrum is used to get the bias b(z) and parameter Q by fitting the auto power spectra of 
the galaxies to the measured values. For the SDSS LRGs, the redshift distribution is determined by the method 
in [561 ] . and the bias is also obtained by adopting the same Q model fitting to the measurements as done for 
the 2MASS samples. For the SDSS quasars, the preliminary redshift distribution is constructed and is used 
to predict the shape of the function /qsoo-i( z ) in the presence of magnification bias. For NVSS, the analytic 
r distribution for fNVSs{ z ) is constrained by NVSS auto power spectrum and cross power spectra with other 
samples. 
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TABLE I: The mean values with la errors and marginalized 95% C. L. of parameter a for the barotropic model using 
different combinations of data sets. 



For each sample, the individual f x (z) need to be recomputed when the cosmological parameters are changed. 
We have modified the code [13] to give the correct background evolution for our case, following Eq. (fH?)) , and 
the second-order differential equation for the total matter density perturbation is also modified due to the 
interaction of the dark matter density perturbation with the vacuum, as given in Eqs. (|48[) or (1531) for our two 
interaction models. The total matter density contrast is given by 

o m {z) = ■ . (61) 

Pdm + Pb 

We have checked that compared with the density perturbation, the dark matter velocity term, 8 dm, on the 
right-hand-side of Eq. (|48|l is small enough to be ignored when we solve the second-order differential equation 
for the barotropic case. 



C. Barotropic model constraints 



We present the constraints for the barotropic model from the combined data sets of CMB+SNIa+BAO, 
CMB+SNIa+LSS(m), CMB+SNIa+glSW, and CMB+SNIa+LSS(b) in TableD Here LSS(b) denotes the con- 
straints obtained by fitting the LSS power spectra to the baryonic matter power spectra, while LSS(m) means 
that the theoretical predictions for the total matter (baryon+dark matter) power spectra are used. The corre- 
sponding ID marginal distribution of the parameter a and the 2D contour plots of the parameters a and Qy 
are shown in Fig. [5] 

From the Table U and Fig. we can see that the joint data from CMB+SNIa+BAO can constrain a to be 
less than or of order 1CP 3 . Tighter bounds on the parameter a can be obtained by using the LSS power spectra 
data. When the total matter power spectra are used to fit the SDSS DR7 power spectra data combined with 
data from CMB+SNIa, the tightest bounds on the parameter a are obtained, to less than or of order 10 -5 . By 
comparison when the baryonic power spectra are are fitted to the SDSS DR7 power spectra data, larger values 
of a are allowed and the upper limit at 95% confidence level (C. L.) is 3.93 x 10 -4 . This is because the dark 
matter power spectra are damped on small scales more quickly than the baryonic power spectra as shown in 
Fig. [7] From the solid lines on the left panel of Fig. [SJ it is seen that the low-multipolc CMB power spectrum 
for the barotropic model is obviously enhanced with the increasing a through the ISW effect. By using the data 
of CMB+SNIa+glSW, we obtain a < 3.97 x 1CT 5 at the 95% C. L.. 



D. Geodesic model constraints 



For the geodesic model, the constraints from the combination of various data are listed in Table [TTJ 
The marginal ID distributions and 2D contour plots from CMB+SNIa+BAO, CMB+SNIa+LSS(m), and 
CMB+SNIa+glSW are shown in Fig. [TO] 

CMB+SNIa+LSS(m) constraints allow negative mean value for a. Compared with results from CMB+SNIa+ 
BAO and CMB+SNIa+LSS(m) constraints, CMB+SNIa+glSW gives a relatively narrow allowed range for 
a. The allowed range of the parameter a at the 95% C. L. in this case are —0.1606 < a < 0.3042 from 
CMB+SNIa+BAO, -0.2111 < a < 0.2069 from CMB+SNIa+LSS(m) and -0.1522 < a < 0.2557 from 
CMB+SNIa+glSW, respectively. By comparison with the results for the barotropic model on the basis of 
the same data sets, the allowed region of parameter a for the geodesic model is obviously enlarged. 

Finally, in Fig. [TT] we show gISW cross-correlation power spectra predicted by using the mean values of 
CMB+SNIa+glSW for the two models. 
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FIG. 9: The ID marginalized distributions for the parameter a and 2D contours with 68% C. L. and 95% C. L. 
for the barotropic model using the combinations of CMB+BAO+SNIa (black line), CMB+SNIa+LSS(m) (red line), 
CMB+SNIa+glSW (blue line) and CMB+SNIa+LSS(b) (magenta line), respectively. The vertical lines on the top 
panels are the 95% confident levels. 
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TABLE II: The mean values with lcr errors for the geodesic model using different combinations of data sets. 



V. CONCLUSIONS 

Any unified dark matter model can be decomposed into dark matter interacting with vacuum energy [30l— f32lj . 
In such a decomposed model, the vacuum energy is coupled to dark matter and the energy-momentum transfer 
between the vacuum energy and the dark matter is given by the gradient of vacuum energy. Different covariant 
forms for the energy-momentum transfer four-vector have different effects on the CMB angular power spectrum 
and LSS power spectrum through the different evolution of inhomogeneous perturbation evolutions, even if 
the background cosmology is the same. Observational constraints on these decomposed models correspond to 
constraints on the interaction between dark matter and the vacuum energy, where we recover ACDM in the 
limit of zero interaction. 

In this paper, we studied the evolution of perturbations in a decomposed Chaplygin gas model and compared 
this with the evolution in a unified Chaplygin gas model with adiabatic perturbations. We considered two 
different models for the energy-momentum transfer in the decomposed case. In one model the energy-momentum 
transfer is along the gradient of the local dark matter density, in which case we recover the evolution equations 
for perturbations of a single barotropic fluid with the same equation of state as the generalised Chaplygin gas. 
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FIG. 10: The ID marginalized distributions for the parameter a and 2D contours with 68% C. L. and 95% C. L. 
for the geodesic model using the combinations of CMB+SNIa+BAO (black line), CMB+SNIa+LSS(m) (red line), and 
CMB+SNIa+glSW (blue line), respectively. The vertical lines on the top panels are the 95% confident levels. 




FIG. 11: The gISW cross-correlation power spectra with different galaxy samples for the barotropic model (blue line) 
with mean values taken from CMB+SNIa+glSW constraints in the Table [T] and the geodesic model (red line) with mean 
values taken from CMB+SNIa+glSW constraints in the Table [TTJ respectively. The black dots with error bars denote 
the gISW data from Ho et al. [371 ]. 
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The pressure perturbations are adiabatic and the rest-frame sound speed equals the adiabatic sound speed. 
This non-zero sound speed has an effect on both the CMB power spectrum and the LSS power spectrum. 
Increasing a enhances the CMB power spectra at the low multipoles. In addition, the effect on the total matter 
(baryon+dark matter) power spectrum caused by the sound speed is that there are decaying oscillations in the 
matter power spectrum for positive a, as for the generalised Chaplygin gas [20| . 

A different possibility which we considered is that the energy-momentum transfer is along the dark matter 
four-velocity, which leads to the dark matter following geodesies. Compared with the barotropic model, the 
geodesic model has a relative entropy perturbation between dark matter and vacuum energy. Thus the sound 
speed of the interacting dark matter+vacuum energy in the dark matter rest frame is zero. There are no 
oscillations or instabilities in the total matter power spectrum for positive or negative a. 

In the decomposed Chaplygin gas, different interaction models can lead to very different evolution of per- 
turbations due to the different sound speeds. We test whether these two kinds of decomposed models can be 
supported or distinguished by current observations. We have constrained the barotropic model and geodesic 
model by using various combinations of data sets, including CMB constraints, type la supernovae, BAO distance 
measurements, large-scale structure and the integrated Sachs- Wolfe effect.. 

For the barotropic model the most stringent constraint on the a parameter is of the order of 10~ 5 from 
the combinations of CMB+SNIa+LSS(m), where the theoretical total matter (baryon+dark matter) power 
spectrum is fitted to the LSS power spectrum data. We conclude that a barotropic model must be extremely 
close to the ACDM model. In contrast the allowed region for the a parameter in the geodesic model is much 
and larger a, and negative a, can also be compatible with the observations. The geodesic model thus allows 
significant deviations from the ACDM model. 

In the case of a unified dark matter model it is not clear whether the galaxy power spectrum on large scales 
should follow the total matter power spectrum, which gives rise to the gravitaitonal potential, or the baryonic 
matter power spectrum, since only baryons follow geodesies. This choice leads to different predictions for the 
LSS power spectrum |22|. The same ambiguity applies in our barotropic model of interacting vacuum+dark 
matter, since the dark matter does not follow geodesies in this case, so we present constraints using LSS fitted to 
both the total matter power spectrum and the baryonic matter power spectrum. In the geodesic model, however, 
we only use the total matter power spectrum to fit the LSS power spectrum since both the dark matter and 
the baryons follow geodesies in this case and we expect the total matter overdensity to be responsible for the 
formation of collapsed halos. Nonetheless it would be interesting to study further the process of non-linear 
collapse in this model and we hope to return to this question in future work. 
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Any unified dark matter cosmology can be decomposed into dark matter interacting with vacuum 
energy, without introducing any additional degrees of freedom. We present observational constraints 
on an interacting vacuum plus dark energy corresponding to a generalised Chaplygin gas cosmology. 
We consider two distinct models for the interaction leading to either a barotropic equation of state 
or dark matter that follows geodesies, corresponding to a rest-frame sound speed equal to the 
adiabatic sound speed or zero sound speed, respectively. For the barotropic model, the most stringent 
constraint on a comes from the combination of CMB+SNIa+LSS(m) gives a < 5.66 x 10 -6 at the 
95% confidence level, which indicates that the barotropic model must be extremely close to the 
ACDM cosmology. For the case where the dark matter follows geodesies, perturbations have zero 
' sound speed, and CMB+SNIa+glSW then gives the much weaker constraint —0.15 < a < 0.26 at 

the 95% confidence level. 

I. INTRODUCTION 

-xT : 

The dominance of dark energy and dark matter in the present cosmological energy density has led many 
/"*\ authors to seek a unified description in terms of a single dark component, referred to as unified dark matter or 
quartessence which should explain both the current accelerated expansion of the Universe and the role of 
nonbaryonic dark matter in structure formation. One candidate from this kind of model, which has attracted a 
' l°t OI attention, is the generalized Chaplygin gas (GCG) model It is a fluid, defined by an exotic equation 

Qh, of state, 

O ' A 

P 9 c a = --^, (1) 

c$ ; 

with model parameters a and A. The GCG model with a = 1 reduces to the original Chaplygin gas model 
' HJ. Predictions of the GCG cosmology has been tested against the different observational data, including 

type la supernovae (SNIa) [5|-tZj], the positions of cosmic microwave background (CMB) radiation peaks [8Hl0j . 
I/-) the X-ray gas mass fraction of clusters [HI , Hubble parameter-redshift data [TH, [l3[ , the lookback time-redshift 
7— I ' data [lH , gamma-ray bursts [l5| , and as well as various combinations of the updated data [l6| . 

Going beyond constraints provided by the homogeneous background cosmology, the inclusion of inhomoge- 
neous linear perturbations in the GCG has an effect on both the cosmic microwave background anisotropies 
and structure formation. In particular, the imprint of a large integrated Sachs- Wolfe (ISW) signal [17j . which 
arises from a time-dependent gravitational potential when dark energy dominates at late times in the Universe, 
could be observed on the CMB anisotropies [l^, [l!| . Recently, in Ref. [2(| two of us revisited constraints on 
the GCG model from the combination of CMB, baryon acoustic oscillations (BAO), and SNIa, including the 
full CMB temperature and polarization power spectra data. A tight bound for the parameter a was obtained 
up to the order of 10~ 3 . 

The effect of the GCG on structure formation also contains important information. In Ref. [2l[, it was argued 
that in order to predict large-scale structure (LSS) consistent with the data from the 2dF survey J22|, the GCG 
model has to be indistinguishable from ACDM, with a < 0(1O~ 5 ). The slightly looser constraint, a < O(10~ 4 ) 
was obtained in Ref. [23| by fitting the predicted baryon matter power spectrum, rather than the total matter 
power spectrum, to the Sloan Digital Sky Survey (SDSS) DR7 power spectrum data 1 . In either case a stringent 
limit is placed on the a parameter because there are either instabilities or oscillations in the GCG power spectra 
due to the adiabatic pressure perturbation produced for nonzero a. This problem can be avoided in some unified 
dark matter models by requiring a fast transition between a CDM-like era and a ACDM- like phase [251427^ 

In order to save the GCG model some type of intrinsic entropy perturbation must be introduced in the GCG 
fluid to counteract the effect caused by the adiabatic GCG pressure perturbation, as in the so-called silent 
GCG model (28|. The combined constraints from data on both the silent GCG and standard GCG model are 
presented in Ref. [29^ . It is shown that with a vanishing GCG pressure perturbation there is a much wider 
allowed parameter range, —0.3 < a < 0.7 [2!|. Subsequently, the ISW data [3(3] was used to constrain the silent 



X 



1 The theoretical baryon power spectrum may also be compatible with SDSS observations for large a > 1, corresponding to a 
superluminal sound speed [23l . however the CMB power spectrum does not favour this regime ^3). 



2 



GCG model 31]. The required entropy perturbation can be introduced by decomposing the GCG fluid into two 
components of dark matter and dark energy [32| , but the possible ways of decomposition are not unique [33| . 
In Ref. |34| the GCG model was split into dark matter and vacuum energy, considering a homogeneous vacuum 
when studying density perturbations in the Newtonian limit. 

Recently it has been shown that any unified dark matter model can be decomposed into dark matter inter- 
acting with a vacuum energy (with equation of state P = —p) without introducing any additional degrees of 
freedom [351 ] . In general an interacting vacuum energy is inhomogenous, and a description of the dynamics of 
inhomogeneous vacuum energy was presented in Refs. [35l. l36j|. 

In this paper we will study in detail the constraints on a decomposed GCG model described by dark matter 
interacting with inhomogeneous vacuum energy. A closed set of perturbation equations are obtained once a 
covariant form is specified for the interaction. Firstly, we consider a model in which the vacuum energy is 
a function of the local matter density. We show this leads to a barotropic equation of state for the dark 
mat ter+ vacuum, and hence reproduces previous results for a unified GCG. Secondly we consider an energy- 
momentum transfer parallel to the matter four-velocity, leading to a geodesic flow for the matter. In both cases 
we consider background cosmologies that mimic the GCG cosmologies, but the evolution of linear perturbations 
differs due to the presence of non-adiabatic pressure perturbations. We pay particular attention to their imprints 
on the CMB and LSS power spectra. To constrain the barotropic and geodesic model parameters we will use 
different combinations of current cosmological data including the CMB temperature and polarization power 
spectra [2Z[, SNIa from the Union2.1 compilation of the Supernova Cosmolo gy P roject Collaboration [38|, the 
BAO distance measurements (39|, the LSS power spectrum from SDSS DR7 [40| . and the galaxy-ISW (gISW) 
cross-correlation power spectra from Ref. [4l| . 

This paper is organized as follows. In the next section, we review the background equations for the GCG 
model. In Sec. Ill, the linear perturbation equations for both the unified GCG model and the two decomposed 
models are presented. In Sec. IV, we show the effects on CMB and LSS power spectra in the two decomposed 
models, then perform a global fitting to observational data by using the Markov Chain Monte Carlo (MCMC) 
method and discuss the resulting constraints on the interaction. In the final section we present our conclusions. 
In what follows (unless otherwise specified) the illustrative plots are shown using the joint mean values for 
cosmological parameters from WMAP7 data [37|: il b h 2 = 0.02255, fl dm h 2 = 0.1126, h = 0.702, r = 0.088, n s = 
0.968, and A s = 2.43 x 10~ 9 at k s0 = 0.002Mpc _1 . 

II. CHAPLYGIN GAS COSMOLOGY 

We will consider a spatially flat, homogeneous and isotropic Friedmann-Robertson- Walker universe with 
metric 

ds 2 = -dt 2 + a 2 (t)Si j dx i dx j , (2) 
where the Friedmann equation is given by 

d\ 2 8nG 



H =(-) = — (Pb + Pr + PgCg) , (3) 

where pb, p r and p g Cg are the energy densities of baryons, radiation and a generalised Chaplygin gas which can 

act as both a dark energy component, accelerating the Universe at late times, and dark matter. The energy 
conservation equations for the various components read 

p b + 3Hp b = , (4) 

p r + AHp r = , (5) 

p gCg + 'SH( PgCg + P gCg ) = . (6) 

The equation of state for GCG {1} then allows us to integrate Eq. © to obtain 



P g c g - P g c g o 



.-il/fl+a) 

B s + (1-B s )a- 3 ^ , (7) 



where p g c g o is the present value of energy density of GCG (when a = do = 1) and we have replaced A in Eq. ([T|) 
with the dimensionless model parameter, B s = A/ 'p g cg*Q- 

In a previous work we have shown that a unified dark matter fluid with an arbitrary equation of state can 
be decomposed into pressureless dark matter interacting with a vacuum energy, V, with an appropriate energy 
transfer [35[ (see also Ref. [II])- For the GCG solution (J7J), we can write 



PgCg — Pdm + V , 



(8) 
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such that, from Eq. ([T]) 

A = V( Pdm + V) a . 
The conservation equations of dark matter and vacuum energy are 



Pdrn + SHpdm — ~ Q ■ 

V = Q. 



Thus from Eqs. ©-(HI]), we obtain 



Q = 3aH 



PdmV 



Pdr, 



V 



(9) 

(10) 
(11) 

(12) 



In the decomposed model (0), the parameter, B s in Eq. (JTJ), can be written in terms of the present values of 
the energy densities, B s — [V/(pd m + V)]o- The decomposed model is characterised by a single dimensionless 
parameter a [35| , but an FRW solution is specified by two boundary conditions (such as the present density of 
matter and vacuum) whereas the original GCG model is specified by two model parameters, a and A (or B s ) 
and one boundary condition. This already suggests that the decomposed model may naturally accommodate 
nonadiabatic density perturbations, in addition to the adiabatic perturbations of the original GCG. As we shall 
see, this can have important consequences for the evolution of perturbations in a decomposed model. 

Note that the Friedmann equation ([3]) can be written as 



H 2 = nl 



6 a 



-3 



+ (Sly + fldm) 



n 



din 



-3(l+a) 



dm 



fly + Qdm 

where the dimensionless fractional energy density of individual components is defined as 



(13) 



(14) 



which satisfies fib + O r + fly + fldm = 1- 

We illustrate the background evolution as a function of redshift for different values of the parameter a in 
Fig. [T] In the particular case when a = 0, the background evolution is identical to that for the ACDM model. 




FIG. 1: The Hubble parameter as a function of redshift z for different values of parameter a, where the solid line with 
a = corresponds to ACDM model. 



The FRW background solution, Eq. (113)) . shows a degeneracy between the decomposed vacuum+dark matter 
with the interaction (fT2"j) and the unified GCG with equation of state ([T]), where £l g cg — + ^dm and 
B a = fly /{Qv + ^dm)- As discussed in Ref. [33j, this is an example of a more general degeneracy between 
interacting and unified models. One can either specify an exotic equation of state for a unified dark matter 
model, or one can obtain many equivalent decomposed models containing interacting fluids with a specified 
energy transfer. Although equivalent to a unified model at the background level, we expect different physical 
models for the interaction to be distinguished by the evolution of inhomogeneous perturbations. 



III. THE LINEAR PERTURBATION EQUATIONS 

We begin with a general scalar mode of the metric perturbation in an arbitrary gauge with the perturbed 
line of element I42H451 



ds 2 = -(1 + 2<j))dt 2 + 2ad l Bdtdx i + a 2 [(1 - 2^)5 ij + 2d l d 3 E) dx i dx j 



(15) 
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The Einstein equations are (42H4 

3H i^J + Hep) -—(ip + Ha) = -4irG6p , (16) 

4> + H(j) = -4irG(p + P)9, (17) 

& + Ha-<p + ip = 8ttGU, (18) 

.... . / 2 V 2 \ 

ip + 3Ht/j + H<t> + {3H 2 + 2H)cf> = AttG ( SP + -—II J , (19) 

where a = a 2 E — aB is the shear, Sp is the total density perturbation, SP is the total pressure perturbation, IT 
is the total anisotropic stress, and 9 is the total covariant velocity perturbation 2 . 

For each individual component, the covariant conservation equation allowing for energy-momentum transfer 
gives 

V^Tft = Qfo . (20) 

General relativity requires conservation of the total energy-momentum, V^T^ = 0, so we have J2i Q(i) = 0- 
The perturbed energy-momentum transfer four- vector is conventionally split into the energy transfer, Qi + SQi, 
and momentum transfer, //, relative to the total four- velocity [44TJ47I ] 

Q { £> = hQ/(l + 4>) - SQiAifi + Qj9)] ■ (21) 
The energy and momentum conservation equations for each fluid component are then [43M45L \at\ : 

5p I + 3H(5p I + 5P I )-3(p I + P I )ij + ( Pl + P I )^r(e I + a) = SQj + Qrf, (22) 

or 

(pi + Pi)6! - 3c 2 sI H( Pl + + (pj + P 7 )0 + SPj + ~m = fi + Qj9 - (1 + c^QjOj , (23) 

3 a z 

where c 2 sI — Pi/ pi is the adiabatic sound speed. 

The numerical code CAMB 48] is written in the synchronous gauge, so in order toplace quantitative con- 
straints on the model parameters we will solve the equations in a synchronous gauge [49J where 

<j) = B = 0, i> = T}, and k 2 E = -- - 3rj . (24) 

There is a residual spatial gauge freedom in the synchronous gauge which allows one to set the initial velocity 
of any component to zero. We see from Eq. (|23p that it is possible to work in a gauge in which 9j — at all 
times for non-interacting, pressureless matter. A comoving synchronous gauge is commonly used, in which the 
dark matter velocity is zero, 9dm = 0. This is no longer possible in models where the dark matter experiences 
a momentum transfer, fd m + Qdm{9 — 9dm) 7^ 0. We can still work in a synchronous gauge, but note that the 
dark matter velocity is not necessarily zero. 

The time-time and time-space components of the perturbed Einstein field equations, (I16|) and (1171) . in the 
synchronous gauge are 

+ (25) 

r) = -AttG(p + P)9. (26) 

Together with the energy and momentum conservation equations for each component, these form a closed set 
of equations. 

After decoupling, the baryon component is conserved independently and the perturbation equations for the 
baryon density contrast and velocity are given by 

k 2 h 
a^° b ~ "2 

9\ = 0. (28) 
Combining Eqs. (fT8"|) . ([T9")) in the synchronous gauge and Eq. (|2l)]) . we have 



h + 2Hh = -8nG(Sp + 3SP) . (29) 

By differentiating Eq. (j2"T)) with respect to time and combining Eq. (|2"8j) and the above equation, we find a 
second-order differential equation for the baryonic density contrast, 

5 b + 2HS b = AnG(6p + 3SP) . (30) 



2 Note that 9 = a(v + B) in terms of the total velocity potential v defined in, e.g., Ref. l4- r l . 
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A. Unified GCG model 



First let us consider perturbations in the unified GCG model, where the energy-momentum tensor of the GCG 
is conserved and hence Eqs. (j2"2")l and ([25]) give the perturbed energy and momentum conservation equations 

( k 2 h\ 

5gCg + 3H(c 2 sgCg -w g cg)5 gCg + (l + wgcg) \--^0 gCg + -I =0, (31) 

r 2 

C s,gCg 



8gCg - 3c c gH9 gC g + — '— Sgcg - , (32) 

1 + W„Co 



w gCg 

where the adiabatic sound speed for the GCG is fixed by the equation of state of the GCG 

c l. q cq = ~ aw gCg, (33) 



Be 

UgCa B s + (1 - B s )a- 3 ( 1 +«) 



(34) 



We can derive a second-order differential equation for GCG overdensity by differentiating Eq. (|3"Tj) with 
respect to time and combining Eqs. (|2T)f and (|3"2"j) . 

SgCg + 2H[l+3(c 2 sgC g-WgCg)}8g C g 

+ [3H(c 2 s gCg - Wgcg) + 3H(c 2 gCg - Wgc g y + m 2 {c 2 s gCg - w gCg ) 2 + &H 2 {c 2 gCg - w gC g)}5gCg 

k 2 k 2 
= - c \gcg^gcg + (1 + w gC g)3clg C gH—0g C g + (1 + w gC g)^G{5p + 36P) . (35) 

The first term on the right-hand-side of the above equation dominates in the small-scale limit and the adiabatic 
pressure perturbation causes the instabilities (or oscillations) of density perturbations for negative (or positive) 
a 501. 



B. Interacting vacuum energy + dark matter models 

Let us now consider a general interacting vacuum energy + dark matter model, where the energy continuity 
equations (|22[) for the dark matter and vacuum energy in the synchronous gauge become 



k 2 „ h 



dPdm + SHSpdm + Pdrn J#dm + ^ I = ^Qdm , (36) 



SV = SQ V , (37) 

while the momentum conservation equations (|23[) are given by 



PdmQdm — fdm + Qdm{d~0dm), (38) 

-SV = fv + Qvd- (39) 

Note that the vacuum energy has vanishing momentum and therefore the force exerted on the vacuum, di(fv + 
Qv&): is exactly canceled by the gradient of the vacuum pressure, SV). 

The conservation of total energy-momentum for the interacting dark matter and vacuum energy implies that 

Qv = —Qdm — Q , SQv — —5Qdm , fv = —fdm ■ (40) 

We can then eliminate SQdm and fdm using the energy and momentum conservation equations for the vacuum 
(|37[) and (|39p to leave the energy-momentum conservation equations for the dark matter 



Q z k 2 Q h SV 

-^Vdm + 77 — 

a 2 p dn 

Pdrn "d m — 

SV + Q0 dm , (42) 



''■In: - Sdm 2® dm o — ~ ~ fety 



where the background energy transfer Q is given by Eq. (fI3 

Note that the interacting vacuum+dark matter introduces no additional degrees of freedom with respect 
to the original GCG. We have two coupled first-order equations (14"TT) and (l4"2"j) for the dark matter density 
and velocity, instead of the two coupled first-order equations (|31|) and (j32|) for the GCG density and velocity. 
The background solution determines the adiabatic sound speed and equation of state for the GCG, whereas a 
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physical model for the interaction is required to determine the vacuum perturbation, 6V, and can affect the 
sound speed for the interacting vacuum+dark matter. 

Generally we can identify two classes of interacting models which could be used to specify the energy- 
momentum transfer in the decomposed model. Either one can specify the vacuum energy as a function, V(X), 
of some scalar quantity X, and hence SV — V'(X)6X. Or else one can identify a four-vector, U^, which the 
energy-transfer (|2"0)) is parallel to: Q v oc C/ M . In the following we will consider one example of each class of 
interacting model, but we note that these are not unique and it may be interesting to consider other possibilities 
with different observational signatures. 



1. Barotropic model 

Firstly, we consider a model in which the vacuum energy density is a function of local dark matter density 

V = V( Pdm ) . (43) 

This implies that the energy transfer between the vacuum and the dark matter is along the gradient of the local 
dark matter density: Q v oc V^pdm- 
At linear order, we have 



sv = 



V 



-Sp 



dm ■ 



(44) 



which shows that in this case the relative perturbations between the vacuum and the dark matter are adiabatic. 
The total pressure perturbation for the interacting dark matter+vacuum energy is 



SV = 



aV 



Pdm + V 
2 



(Spdm + SV) 



gCg( 6 Pdm + SV) . 



(45) 



Therefore, the sound speed for the interacting dark matter+vacuum energy is exactl y th e same as the adiabatic 
sound speed for the original generalised Chaplygin gas, c 2 gCg given by Eq. (j3"31 [35|,|36j]. We show the effect on 
the sound speed caused by the parameter a in Fig. [2] 



FIG. 2: The sound speed as a function of scale factor a for the barotropic model. 



In the synchronous gauge, we have from Eqs. (|4"T]) . (l4"2l and (|4"4"|) 



™drn 



(l + a)Q p dm + (l + a)Vk 2 

-Odrn H — Tr o^dr, 



Pdn 



Pdm + (1 + a)V 



V 



Pdn 



V 



a)V h 
2 ' 



_Q_ 

Pdm 



aV 



dvn 



p dm + (1 + a)V 



'dm ■ 



(46) 
(47) 



It is clear that the barotropic model with a = reduces to the ACDM model at the linear order. Comparing 
with dark matter perturbations in the standard ACDM model, we can see that density contrasts and velocities 
of the interacting dark matter evolve very differently. In this case, due to the presence of momentum transfer, 
Qdm will evolve to be nonzero, even if its initial value is set to zero. Therefore, the synchronous gauge is not 
comoving with the interacting dark matter and we calculate the perturbed equations in an arbitrary synchronous 
gauge [51| . 
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Combining Eqs. (g5) and 07} with Eq. 
matter overdensity 



we obtain a second-order differential equation for the dark 



-din 



2(1 + a)Q 



Pdm + (1 + 
(l + a)Q 



2H 



^dm 



(l + a)Q 



p dm + (1 + a)V\ [p dm + (1 + a)V 
k\ f i | «V \ Q fc 2 f 

\ Pdm + V / Pdm a 2 



2H 



(l + a)Q 



C s,gCg n 2^ dm 



1 



Pdm + (1 + a)V 



3dm 



Pdm + ^ 



47rG((5p + 35P) . 



(48) 



The first term on the right-hand-side of Eq. (|48|) dominates in the small-scale limit. Negative a causes insta- 
bilities of density perturbations. Conversely, there are oscillations of perturbations for positive a. In order to 
avoid the presence of negative sound speed, we will require a > in this case. 

We find that the perturbations for the barotropic model, Eqs. (I4l)l) - (j4"5|) . obey the same dynamical evolutions 
as those for the unified GCG model, Eqs. (|31[) . (l32l and (13"5j) . where we identify 



S gCg - 



Pdn 



+ {l + a)V 



y gCg 



J dm • 



(49) 



We should not be surprised that this interacting model reproduces the same behaviour as the GCG since by 
requiring V — V{pdm) we have set the pressure of the interacting mattcr+vacuum P dm +v = —V to be an 
implicit function of pdm + V, i.e., we have a barotropic fluid. 

The evolution of the dark matter density perturbations for the barotropic model on a fixed scale is shown in 
Fig. [3] In the left panel, it is shown that increasing a causes the suppression of the growth of density fluctuations 
on the scale, k — 0.001 [h Mpc -1 ]. For the smaller scale, k = 0.01 [h Mpc -1 ], in the right panel, we can see that 
density fluctuations rapidly decay and oscillate with increasing a. 





0.2 0.4 0.6 



FIG. 3: The dark matter density perturbations on fixed scales: k = 0.001 [h Mpc -1 ] (left panel) and k = 0.01 [h Mpc -1 ] 
(right panel) as a function of scale factor for the barotropic model. The thin solid line, dashed line, dotted line, dot-dashed 
line and thick solid line correspond to a — 0.1, 0.05, 0.01, 0.001, (ACDM model), respectively. 



2. Geodesic model 



An alternative model for the energy transfer is to consider an energy flow that is parallel to the four- velocity 
of the dark matter: 



-(dm) — Q U (dm 



(50) 



For this form of the covariant interaction there is no momentum transfer in the rest frame of dark matter 
[13) HI]) fdm — Q(8 — 9dm) = in Eq. (p?T]) . The velocity perturbation for dark matter is not affected by the 
interaction and obeys the standard equation, 9 dm = from Eq. (|38p . This allows us to work in a synchronous 
gauge that is comoving with the dark matter (0dm = 0). The dark matter four- velocity, is thus a geodesic 
flow: u^dm)^ ^(dm) = ^- Hence we will refer to this form of energy-transfer as a geodesic model. 
The vacuum perturbation in the dark matter-comoving frame is identically zero, from Eq. (|42|) . 



SV com = SV + V9 dm = . (51) 
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The pressure perturbation for the interacting vacuum+dark matter is therefore zero in the dark matter rest 
frame and thus the speed of sound for the interacting vacuum+dark matter is zero 3 [III HH . Note that this 
implies that in this case there is a non-adiabatic pressure perturbation, due to the relative entropy perturbation 
between the dark matter and the vacuum energy (35l . [3(| , in contrast to the barotropic model. 

From Eq. (|4ip , the perturbation equation for the dark matter density contrast in the rest frame of dark matter 
obeys 



5dm — —- 



Pdn 



-&dn 



(52) 



We see that, again, the case with a = is completely equivalent to ACDM at linear order. 

A second-order differential equation for the dark matter overdensity can be derived from the above equation 
and Eq. (j2"!)|) , which gives 



J dm 



--+2H )<),,„ 

Pdm 



2H 



Q_ 

Pdm 



inG(5p + 35P) . 



(53) 



We plot the evolution of the dark matter density perturbations on different scales in Fig. 01 Compared 
with the growth in the ACDM model, positive a indicates that there is an energy transfer from dark matter 
to vacuum energy, so the growth of density perturbations with positive a is suppressed at early times during 
matter-domination. At late times, vacuum energy dominates in the Universe, and the vacuum interaction 
drives the growth of dark matter density perturbations. Conversely, for the negative a the growth of density 
perturbations is enhanced at early times and is gradually suppressed at late times. Note that there is no 
instability associated with imaginary speed of sound for a < in the geodesic model. 






-.0.2 


(X= 


=0.1 




=0 
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FIG. 4: The dark matter density perturbations on fixed scales: k = 0.001 [h Mpc -1 ] (left panel) and k = 0.01 [h Mpc -1 ] 
(right panel) as a function of scale factor for the geodesic model. The thin solid line, dashed line, thick solid line, dotted 
line and dot-dashed line correspond to a — 0.2, 0.1, (ACDM model), -0.1, -0.2, respectively. 



IV. COSMOLOGICAL CONSTRAINTS 



A. Theoretical predictions of CMB and LSS power spectra 

First, we illustrate the imprint on the CMB and LSS power spectra for different values of the parameter a 
for both the barotropic model and geodesic model. For the barotropic model we will restrict the value of a to 
be greater than or equal to zero in order to avoid instabilities due to an imaginary sound speed, but we will 
extend a to negative values for the geodesic model since the sound speed remains real in that case. 

The effects on the CMB power spectra are shown in Fig. [5] From the left panel, positive a leads to an 
energy transfer from dark matter to vacuum energy. Increasing a causes the matter-radiation equality to 



A vanishing pressure perturbation was also considered in a Cardassian model without any internal fluctuations, which in a 
special case that v = 1 reproduces a GCG background 54|. However we will show that in the geodesic model dark matter density 
perturbations are affected by the background interaction and obey a modified perturbation equation. 
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happen earlier and hence suppresses the first acoustic peaks. Conversely, negative a delays the epoch of matter- 
radiation equality, and hence amplifies the acoustic peaks, as shown in the geodesic case in the right panel of 
Fig. [5] 

One source of CMB anisotropics at low multipoles is the ISW effect. The amplitude of the ISW signal on large 
scales is affected not only by the epoch of dark matter-vacuum energy equality, but is also modified directly 
by the evolution of the gravitational potential, which differs in the barotropic and geodesic models. These two 
competing mechanisms give the final ISW signal. 

In order to clearly express the ISW effect caused by the time- varying gravitational potential, we calculate the 
Newtonian potential, i/j , corresponding to the metric perturbation, tp, in the Newtonian or longitudinal gauge 

(B = E = 0) m. 

This is given by the relativistic Poisson equation [44| 

k 2 ^ N = -47rGa 2 pA , (54) 
where the total comoving density perturbation is given by 

pA = J2p I [5 I -3H(l + wi)ei] , (55) 

including the vacuum perturbation. The traceless, space-space component of the Einstein field equations, (fT8l) . 
gives the relation between the two scalar potentials in the Newtonian gauge [44| : 

■4> N - <j> N = 8irGH , (56) 

where II is related to total the anisotropic stresses. 

Now the ISW temperature anisotropy is given by the following source: 



N i IN 



Sisw = ip 

= - 8 ^(^ A + n )" < 57 > 

At late times, the relativistic particles (photons and neutrinos) can be neglected and there is no anisotropic 
stress. Therefore, it is the comoving density perturbation that drives the late-time ISW effect. Figure [5] shows 
how the parameter a affects the time evolution of the gravitational potential in the two models. We see that for 
the same value of a we find a much larger variation in the Newtonian potential at late times in the barotropic 
model due to the nonzero sound speed, as shown in Fig. [5] This will be reflected in much tighter bounds on the 
value of a in the barotropic model coming from the overall CMB anisotropics, and these constraints become 
even tighter when the galaxy clustering data is also used to identify the late-time ISW effect. 

Note that the overall CMB anisotropics, including the ISW effect, are the same in the barotropic model and 
the original unified Chaplygin gas model 4 since the evolution of the total matter+vacuum density perturbations 
is the same as the density perturbations in the unified GCG model. 




FIG. 5: CMB temperature power spectra vs multipole moment 1. The solid line is for the ACDM model. Left panel: The 
thick lines and thin lines for the same values of parameter a correspond to the barotropic model and geodesic model, 
respectively. Right panel: The dashed line and dotted line correspond to a = 0.1 and a = —0.1 in the geodesic model, 
respectively. 



4 This case can also be obtained from case II of the Cardassian model in Ref. [54| for v = 1. 
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FIG. 6: The evolution of the gravitational potential as a function of scale factor on the large scale k = 0.001[h Mpc -1 ]. 
The thick solid line is for the ACDM model. The thick dashed, dotted and dot-dashed lines correspond to a = 0.1, 0.2, 
0.3 for the barotropic model, respectively. The thin dashed, dotted and dot-dashed lines correspond to a = 0.1, 0.2, 0.3 
for the geodesic model, respectively. The thin solid line is the evolution of the gravitational potential for a — —0.1 in 
the geodesic model. 



In Fig. [7J we present the dark matter power spectrum, the baryonic power spectrum, and the total matter 
power spectrum (baryon and dark matter) for the barotropic model. The figure shows that, as expected, the 
dark matter has an oscillating power spectrum in the small-scale limit for nonzero a. The larger the parameter 
a, the larger the scale where the oscillations of dark matter power spectra begin. By comparison, the baryonic 
power spectra behave much more smoothly as they are only gravitationally coupled to the oscillating dark 
matter density (23j . The total matter power spectra combining both dark matter and baryon densities therefore 
shows intermediate behaviour, but cannot avoid the rapid decay of power spectrum on the small scales for 
nonzero a. 

Figure [8] shows the total matter power spectra for the geodesic model. The problematic oscillations of the 
dark matter power spectra do not occur in this case since the sound speed is zero. The only oscillations seen are 
the baryon acoustic oscillations, also seen in ACDM. The matter power spectrum is enhanced on small scales 
for a > due to the earlier matter-radiation equality which moves the turnover in the matter power spectrum 
to smaller scales. 



B. Constraint method and data 



Next, we will use a variety of cosmological data to constrain these two models for the matter+ vacuum 
interaction, according to the Markov Chain Monte Carlo method, using the publicly available CosmoMC package 
55, 56). The cosmological parameter set is P = {flbh 2 , SI,/,,,/; 2 , 0g, t, o . ??,,, log[10 10 A s ]}, where Sl b h 2 and fl,/ m /i 2 
are the physical baryon and dark matter densities, 0s is the ratio (multiplied by 100) of the sound horizon 
and angular diameter distance, r is the optical depth, a is the dimensionless generalised Chaplygin gas model 
parameter, n s is the scalar spectral index, and A s is defined as the amplitude of the initial power spectrum. 
Note that the dimensional parameter A in the original definition of the Chaplygin gas ([I]) is given in terms of 
fldm and fly using Eq. §§§ as 

A = n v (n dm + n v ) a , (58) 

and fly = 1 — Ob — 57 r — Qdmi by the requirement of spatial flatness. 

We will use different combinations of data (CMB+SNIa+BAO, CMB+SNIa+LSS, and CMB+SNIa+glSW). 
For the LSS power spectra we will only use the linear matter power spectrum up to k ~ 0.1[h Mpc -1 ]. In 
addition, we use a top- hat prior for the cosmic age, i.e., lOGyr < to < 20Gyr, impose a weak Gaussian prior on 
the physical baryon density fi{,/i 2 = 0.022 ± 0.002 from big bang nucleosynthesis [13], and use the value of the 
Hubble constant Hq = 73.8 ± 2.4km s _1 Mpc _1 [58[ taking a Gaussian likelihood function. 

The theoretical galaxy-ISW (gISW) cross-correlated power spectrum, P s $, where $ = ip N + cb N , can be 
calculated by a modified CAMB code. We modify the code by referring to ISiTGR [59l| . so that a power 
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FIG. 7: The LSS power spectra as a function of wavenumber for the barotropic model. In each panel, the solid line is the 
total matter power spectrum for the ACDM model. The dashed line, dotted line and dot-dashed line are the baryonic 
power spectrum, the total matter power spectrum and the dark matter power spectrum, respectively. The separate 
panels (left to right and top to bottom) correspond to a — 10 -5 , 10 -4 , 10 -3 , 10 -2 , respectively. 
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FIG. 8: The total matter power spectra as a function of wavenumber for the geodesic model. The thin solid line, dashed 
line, thick solid line, dotted line and dot-dashed line correspond to a = 0.2, 0.1, (ACDM model), -0.1, -0.2, respectively. 



spectrum can be directly calculated from the individual transfer function. In the absence of magnification bias, 
the gISW cross correlation is written as 



rSx T = T C mb 
1 (/ + 1/2) 2 



(/:/),.( : Ill, ( : )P. A , ( - ^*/ 2 , z 



(59) 



where the galaxy samples in Ref. [4JJ are collected from the Two Micron All Sky Survey (2MASS) Extended 
Source Catalog (XSC), which is divided into 4 bins: 2MASS0-3, the luminous red galaxies (LRGs) from SDSS 
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split into 2 sample bins: LRGO-1, the photometric quasars from SDSS with 2 redshift bins: QSOO-1, and NRAO 
VLA Sky Survey (NVSS): NVSS. Namely, there are nine galaxy sample bins x. The scale-independent bias 
factor b x (z) is assumed to relate the observed projected galaxy overdensity to the matter density. H x (z) is the 
redshift distribution. \( z ) is the comoving radial distance. In ISWWLL code, the function f x (z) = b x (z)Tl x (z) 
is determined by fitting auto power spectra and cross power spectra between the samples: 

Cp* = / dz.f x (z)f y (z)^D\z)P s , (^p.o) , (60) 

even if the methods differ for the various samples. 

For the 2MASS galaxies, the redshifts are identified by matching 2MASS galaxies with SDSS main galaxies. 
The nonlinear Q model [6(| (denoted by the bold text to be distinguished from the background interaction Q) 
of the galaxy power spectrum is used to get the bias b(z) and parameter Q by fitting the auto power spectra of 
the galaxies to the measured values. For the SDSS LRGs, the redshift distribution is determined by the method 
in Ref. [61(, and the bias is also obtained by adopting the same Q model fitting to the measurements as done 
for the 2MASS samples. For the SDSS quasars, the preliminary redshift distribution is constructed and is used 
to predict the shape of the function fQsoo-i(z) m the presence of magnification bias. For NVSS, the analytic 
r distribution for fNVSs( z ) is constrained by NVSS auto power spectrum and cross power spectra with other 
samples. 

For each sample, the individual f x (z) need to be recomputed when the cosmological parameters are changed. 
We have modified the code [4l[ to give the correct background evolution for our case, following Eq. (fT3")l . and 
the second-order differential equation for the total matter density perturbation is also modified due to the 
interaction of the dark matter density perturbation with the vacuum, as given in Eqs. ([35)1 or (15"3l for our two 
interaction models. The total matter density contrast is given by 

r i n Pdmddm + PbSb 

o m (z) = ; . (61) 

Pdm + Pb 

We have checked that compared with the density perturbation, the dark matter velocity term, 9 dm, on the 
right-hand-side of Eq. (135)1 is small enough to be ignored when we solve the second-order differential equation 
(|48p for the barotropic case. 

C. Barotropic model constraints 

We present the constraints for the barotropic model from the combined data sets of CMB+SNIa+BAO, 
CMB+SNIa+LSS(m), CMB+SNIa+glSW, and CMB+SNIa+LSS(b) in Table UJ Here LSS(b) denotes the con- 
straints obtained by fitting the LSS power spectra to the baryonic matter power spectra, while LSS(m) means 
that the theoretical predictions for the total matter (baryon+dark matter) power spectra are used. The corre- 
sponding ID marginal distribution of the parameter a and the 2D contour plots of the parameters a and Qy 
are shown in Fig. [5] 

From Table [J and Fig. [5) we can see that the joint data from CMB+SNIa+BAO can constrain a to be less 
than or of order 1CV 3 . Tighter bounds on the parameter a can be obtained by using the LSS power spectra 
data. When the total matter power spectra are used to fit the SDSS DR7 power spectra data combined with 
data from CMB+SNIa, the tightest bounds on the parameter a are obtained, to less than or of order 10 -6 . By 
comparison when the baryonic power spectra are are fitted to the SDSS DR7 power spectra data, larger values 
of a are allowed and the upper limit at 95% C. L. is 1.07 x 1CV 4 . This is because the dark matter power spectra 
are damped on small scales more quickly than the baryonic power spectra as shown in Fig. [7) From the thick 
lines in the left panel of Fig. [5j it is seen that the low-multipole CMB power spectrum for the barotropic model 
is obviously enhanced with the increasing a through the ISW effect. By using the data of CMB+SNIa+glSW, 
we obtain a < 3.97 x 10" 5 at the 95% C. L. 



D. Geodesic model constraints 

For the geodesic model, the constraints from the combination of various data are listed in Table [TTJ 
The marginal ID distributions and 2D contour plots from CMB+SNIa+BAO, CMB+SNIa+LSS(m), and 
CMB+SNIa+glSW are shown in Fig. [TO] 

CMB+SNIa+LSS(m) constraints allow negative mean value for a. Compared with results from 
CMB+SNIa+BAO and CMB+SNIa+LSS(m) constraints, CMB+SNIa+glSW gives a relatively narrow allowed 
range for a. The allowed range of the parameter a at the 95% C. L. in this case are —0.16 < a < 0.30 
from CMB+SNIa+BAO, -0.20 < a < 0.21 from CMB+SNIa+LSS(m) and -0.15 < a < 0.26 from 
CMB+SNIa+glSW, respectively. By comparison with the results for the barotropic model on the basis of 
the same data sets, the allowed region of parameter a for the geodesic model is obviously enlarged. 
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TABLE I: The mean values with la errors and marginalized 95% C. L. of parameter a for the barotropic model using 
different combinations of data sets. 




FIG. 9: The ID marginalized distributions for the parameter a and 2D contours with 68% C. L. and 95% C. L. for the 
barotropic model using the combinations of data (left to right): CMB+BAO+SNIa (black line), CMB+SNIa+LSS(m) 
(red line), CMB+SNIa+glSW (blue line) and CMB+SNIa+LSS(b) (magenta line), respectively. The vertical lines in 
the top panels are the 95% confident levels. 



Finally, in Fig. [TT] we show gISW cross-correlation power spectra predicted by using the mean values of 
CMB+SNIa+glSW for the two models. 



V. CONCLUSIONS 

Any unified dark matter model can be decomposed into dark matter interacting with vacuum energy (33-l36| . 
In such a decomposed model, the vacuum energy is coupled to dark matter and the energy-momentum transfer 
between the vacuum energy and the dark matter is given by the gradient of vacuum energy. Different covariant 
forms for the energy-momentum transfer four-vector have different effects on the CMB angular power spectrum 
and LSS power spectrum through the different evolution of inhomogeneous perturbation evolutions, even if 
the background cosmology is the same. Observational constraints on these decomposed models correspond to 
constraints on the interaction between dark matter and the vacuum energy, where we recover ACDM in the 
limit of zero interaction. 

In this paper, we studied the evolution of perturbations in a decomposed Chaplygin gas model and compared 
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FIG. 10: The ID marginalized distributions for the parameter a and 2D contours with 68% C. L. and 95% C. L. for 
the geodesic model using the combinations of data (left to right): CMB+SNIa+BAO (black line), CMB+SNIa+LSS(m) 
(red line), and CMB+SNIa+glSW (blue line), respectively. The vertical lines in the top panels are the 95% confident 
levels. 
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FIG. 11: The gISW cross-correlation power spectra with different galaxy samples for the barotropic model (dashed 
line) with mean values taken from CMB+SNIa+glSW constraints in the Tableland the geodesic model (solid line) with 
mean values taken from CMB+SNIa+glSW constraints in the Table [ITJ respectively. The dots with error bars denote 
the gISW data from Ho et al. HJ]. 
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TABLE II: The mean values with la errors for the geodesic model using different combinations of data sets. 



this with the evolution in a unified Chaplygin gas model with adiabatic perturbations. We considered two 
different models for the energy-momentum transfer in the decomposed case. In one model the energy-momentum 
transfer is along the gradient of the local dark matter density, in which case we recover the evolution equations 
for perturbations of a single barotropic fluid with the same equation of state as the generalised Chaplygin gas. 
The pressure perturbations are adiabatic and the rest-frame sound speed equals the adiabatic sound speed. This 
nonzero sound speed has an effect on both the CMB power spectrum and the LSS power spectrum. Increasing a 
enhances the CMB power spectra at the low multipoles. In addition, the effect on the total matter (baryon+dark 
matter) power spectrum caused by the sound speed is that there are decaying oscillations in the matter power 
spectrum for positive a, as for the generalised Chaplygin gas [211 ] . 

A different possibility which we considered is that the energy-momentum transfer is along the dark matter 
four-velocity, which leads to the dark matter following geodesies. Compared with the barotropic model, the 
geodesic model has a relative entropy perturbation between dark matter and vacuum energy. Thus the sound 
speed of the interacting dark matter+vacuum energy in the dark matter rest frame is zero. There are no 
oscillations or instabilities in the total matter power spectrum for positive or negative a. 

In the decomposed Chaplygin gas, different interaction models can lead to very different evolutions of per- 
turbations due to the different sound speeds. We test whether these two kinds of decomposed models can be 
supported or distinguished by current observations. We have constrained the barotropic model and geodesic 
model by using various combinations of data sets, including CMB constraints, type la supernovae, BAO distance 
measurements, large-scale structure and the integrated Sachs- Wolfe effect. 

For the barotropic model the most stringent constraint on the a parameter is of the order of 10 -6 from 
the combinations of CMB+SNIa+LSS(m), where the theoretical total matter (baryon+dark matter) power 
spectrum is fitted to the LSS power spectrum data. We conclude that a barotropic model must be extremely 
close to the ACDM model. In contrast the allowed region for the a parameter in the geodesic model is much 
larger, and negative a can also be compatible with the observations. The geodesic model thus allows significant 
deviations from the ACDM model. 

In the case of a unified dark matter model it is not clear whether the galaxy power spectrum on large scales 
should follow the total matter power spectrum, which gives rise to the gravitaitonal potential, or the baryonic 
matter power spectrum, since only baryons follow geodesies. This choice leads to different predictions for the 
LSS power spectrum [23|]. The same ambiguity applies in our barotropic model of interacting vacuum+dark 
matter, since the dark matter does not follow geodesies in this case, so we present constraints using LSS fitted to 
both the total matter power spectrum and the baryonic matter power spectrum. In the geodesic model, however, 
we only use the total matter power spectrum to fit the LSS power spectrum since both the dark matter and 
the baryons follow geodesies in this case and we expect the total matter overdensity to be responsible for the 
formation of collapsed halos. Nonetheless it would be interesting to study further the process of nonlinear 
collapse in this model and we hope to return to this question in future work. 
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